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In this paper we give a review of the method of imsets introduced by Studenjffl 
from a geometric point of view. Elementary imsets span a polyhedral cone and 
its dual cone is the cone of supermodular functions. We review basic facts on 
the structure of these cones. Then we derive some new results on the following 
topics: i) extreme rays of the cone of standardized supermodular functions, ii) 
faces of the cones, iii) small relations among elementary imsets, and iv) some 
computational results on Markov basis for the toric ideal defined by elementary 
imsets. 



1. Introduction 

The method of imsets by Studeny^ provides a very powerful algebraic 
method for describing conditional independence relations under a proba- 
bility measure. Rules for deriving conditional independence relations are 
translated into relations among integer vectors called imsets. Hence many 
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properties of conditional independence relations can be conveniently in- 
terpreted from a geometric viewpoint. In recent papers Studeny and his 
collaborators^! further develop geometric methods for learning Bayesian 
networks. In this paper we are more concerned on basic geometric proper- 
ties of imsets, in particular from the viewpoint of lattice bases and Markov 
bases for the configuration of elementary imsets. 

The cone of supermodular functions, which we call the supermodular 
cone, is defined by a set of linear inequalities and effective inequalities corre- 
spond to elementary imsets. Hence the ii-representation (cf. Griinbaum 4 ) 
of the supermodular cone is explicitly given. The cone generated by the 
elementary imsets, which we call the imset cone, is the dual to the super- 
modular cone and its set of extreme rays is given by elementary imsets. 
Hence, in the dual sense, the ^-representation of the imset cone is given. 

From an algorithmic viewpoint of convex geometry, the V- 
representation of the supermodular cone, or equivalently the H- 
representation of the imset cone, are hard to compute and characterize. 
Therefore general results on the facets of the imset cone, or equivalently 
the "extreme rays" of the supermodular cone are important. Since the su- 
permodular cone contains a linear subspace consisting of modular functions, 
we consider extreme rays of the cone of standardized supermodular func- 
tions (see ([6]) below for standardization). These extreme rays are called 
skeletal supermodular functions in StudenyP^ 

Although the complete description of the supermodular cone and the 
imset cone is very difficult, some faces of these cones can be studied in 
detail. In particular a face corresponding to a semi-elementary imset seems 
to have a simpler structure than other faces. 

The organization of the paper is as follows. In Sections [2] and [3] we 
present a review of the method of imsets from a geometric point of view. In 
Section [2] we summarize basic facts on the cones of supermodular functions 
and semi-elementary imsets. In Section [3] we review derivations of condi- 
tional independence statements based on multiinformation and imsets. In 
Section 0] we give some results on extreme rays of the cone of standardized 
supermodular functions. In Section [S] we discuss properties of a face corre- 
sponding to a semi- elementary imset. In Section|5]we characterize relations 
among a small number of elementary imsets. Finally in Scction[7]we discuss 
some computational results on Markov basis for the toric ideal defined by 
the configuration of elementary imsets. 



3 



2. Basic facts on supermodular functions and imsets 

As in Kashimura and TakemuraP, we use the notation and definitions from 
StudenyP Let TV be a finite set and let V{N) = {A : A C N} denote its 
power set. In this paper A C B means that A is a subset of B and A C B 
means that A is a proper subset of B. \A\ denotes the cardinality of A. For 
notational convenience, we write the union A U B as AB. A singleton set 
{i} is simply written as i. M., R + , Q, Q + , Z, Z + , N, denote the sets of reals, 
non-negative reals, rationals, non-negative rationals, integers, non-negative 
integers and positive integers, respectively. 
/ : V(N) — > K is called supermodular if 

f(EF) + f(E D F) > f(E) + f(F), VE,FCN. (1) 

The set of supermodular functions over N is denoted by fC(N). f is sub- 
modular if — / is supermodular. 

/ is modular if it is both supermodular and submodular, i.e. 

f(EF) + f(E f) F) = f(E) + f(F), VE,FCN. (2) 

£(N) = JC(N) n (—JC(N)) denotes the set of modular functions over N. A 
modular function / is like a discrete (signed) measure. Indeed if /(0) = 0, 
then by taking disjoint E and F in J5J), we see that / is a measure and 
hence / can be written as f(E) = J2 e eE f( e )- Without the restriction of 
/(0) = 0, we have f(E) - /(0) = E ee£ (/(e) - /W), which can also be 
written as 

f(E) = \$ + J2 A - A » = /W> X - = /( e ) - / W> e e ^ ( 3 ) 

This shows that the dimension of the linear space C(N) is |JV| + 1 and a 
basis of £(iV) is given by the following |JV| + 1 functions: 

h(E) = l, f e (E) = l eE .(E) = l e&E , eeN, (4) 

where l eem (E) = l{ e }c»(E) = l e< zE is the indicator function 

I otherwise. 

Given a supermodular function /, define a modular function and a 
supermodular function / by 

/ L (i?) = A + ^A e , A =/(0), A e = /(e) - /(0), 

eS-E 

f(E) = f(E)-f L (E), (6) 
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Then / is supermodular and f(E) — for \E\ < 1. / is often a preferred 
standardization of /. By induction on the cardinality of \E\, it is easy to 
show that any standardized supermodular function / is non-negative and 
non-decreasing: 

f(E) > 0, ECF^ f(E) < f(F), E, F C N. 

Let N = {1, ... , \N\}. We can identify / : V(N) K with a vector in 
^\-P{N)\ = M 2i"i by listing its va i U e S : 

(/(0),/(l),...,/(|7V|),/({l,2}),...,/(7V)). 

With this identification the set of supermodular functions fC(N) is a poly- 
hedral cone in R 2 '™ , because it is defined by the set of linear inequalities 
in ©. By the definition of £{N) = JC(N) n (-K(N)), £(N) is the largest 
linear subspace contained in K(N). 

For A C N, 5 A ■■ V(N) -> M is defined as 



<Ja(5) 



1 if 5 = A, 
otherwise. 



For pairwise disjoint subsets, A, £?, C C AT, we write this triplet by 
(A, B\C), and the set of all disjoint triplets (A, B | C) over N by T(iV). 
Unless otherwise stated, we assume that A, B are non-empty. On the other 
hand C may well be an empty set. For a triplet (A, B | C) 6 T(N), a semi- 
elementary imset u^a.b | c) : 'P(N) — ¥ K is defined as 

u <a,s | c) = & abc + Sc — Sac ~ $bc- (7) 

If A = or = 0, then uia,b \ c) is the zer0 imset. We consider u^a.b \ c) 
as a 2' Ar l-dimensional integer vector with two elements equal to I (at ABC 
and C) and two elements equal to —I (at AC and BC). If A — a and 
B = b are singletons, the imset w/ a ,6| c) is called elementary. The set of all 
elementary imsets is denoted by £{N). Note that the number of elementary 
imsets is given by 

\£{N)\ - ^h 2 W- a - |JV|(|JV| - l)2l jy l- 3 . 

It is instructive to write out all elementary imsets in a 2^1 x |£(iV)| 
matrix Un, where each elementary imset is a column of Un- For \N\ = 4, 
N = abed, the matrix Un is written in Table [T] In this paper we consider 
this matrix as the configuration defining a toric ideal. Let f(S) = ^\S\ 2 , 
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Table 1. Configuration of elementary imsets Upf, \N\ =4. 
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S C jV. It is easily seen that this / is supermodular and 

J2f(S)u(S) = l, Vue£(N). (8) 

SCAT 

Hence Un defines a homogeneous toric ideal (Chapter 4 of Sturmfels^) . 
Note that since l^l is modular, any f(S) — ^\S\ 2 + ci\S\ + C2, ci,C2 € R, 
satisfies (|8|). In particular f(S) — i\S\(\S\ — 1) is the standard supermodular 
function satisfying ([8]). 

In Table Q] rows are ordered (from bottom to top) by the cardinality of 
the set and by reverse lexicographic order among sets of the same cardinal- 
ity. We call this order the graded reverse lexicographic order of V(N). This 
order can be generalized to elementary imsets as follows. In an elementary 
imset U( a .b | c) we always order a, b as a < b. Then we define 

{C < C, 
or C = & and b < b', (9) 
or C ~ C and b — b' and a < a'. 

In particular, strict inclusion CcC" implies C < C. For N = {1, . . . ,\N\}, 
the maximum element of £(N) with respect to < is un t 2\m, where R = 
{3,...,|iV|}. We call this order the graded reverse lexicographic order of 
£(N). This order will be used in the proof of Theorem 16. II in Section [6l The 
columns in Table Q] are ordered (from right to left) according to the graded 
reverse lexicographic order for elementary imsets. 
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The following identity among semi- elementary imsets is of basic impor- 
tance. 

U(A,BD | C) = V>(A,B | C) + U>(A,D \ BC) ■ (10) 

This identity can be directly verified as 

U (A,B | C) + U (A,D | BC) — ($ABC + $C — &AC — &Bc) 

+ {&ABCD + Sbc — Sabc — 8bcd) 
= Sabcd + Sc — 5 AC — &BCD- 
It is conveniently depicted in Fig. [TJ 

ABCD 




Fig. 1. Sum of two semi-elementary imsets 

From (|10[) we can split A or B in u/a,b\c) into smaller subsets. If we 
repeat this splitting, every semi-elementary imset can be written as a non- 
negative integer combination of elementary imsets 

u (am\c) = X! kv ' v > where fc„ e Z + . (11) 

ves(N) 

In Kashimura et alP we gave a detailed study of the set of all possible 
non-negative integer combinations of elementary imsets which are equal to 
a semi- elementary imset. 

We consider R 2 as equipped with the standard inner product (•,•)■ 
Then the inner product of / : V{N) — > R 2 ' ' and uia,b \ c) lS written as 

(/, u {a ,b | c)) = f(ABC) - f(C) + f(AC) + f(BC). 

This inner product was already considered in (j8]). Let E — AC, F ~ BC, En 
F = C, EF = ABC in (JTJ) . We see that / is supermodular if and only if 

(/,«<A,B|t7>)>0 ) V(A,B\C)eT(N). (12) 
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Hence the set of semi-elementary imsets {u^a,b\c) I B \ C) G T(N)} 
give the ^-representation of the supermodular cone IC(N). By definition, 
the convex cone 

IC*(N) = R+{u {a ,b | c) I (A, B\C)G T(N)} 

= { 51 k( A ,B | C) ' U (A,B\ C) I k( A ,B \ C) > 0} 

| C)dT(N) 

generated by the semi-elementary imsets is the cone dual to the supermod- 
ular cone K.{N). We call K.*(N) the imset cone. 

In the iJ-representation of IC(N), not all of the hyperplanes determined 
by U(a,b\c) are effective. In fact by (fTT|) . if (/, v) > for all elementary 
v € £(N), then (/, U(a,b | c)) ^ for every semi-elementary imset. Hence 
/ is supermodular if (/, v) > 0, W S It also follows that JC*(N) is 

generated by the elementary imsets: 

/C*(iV) =M+£(A0. (13) 

A^* ( A^) is a pointed cone by (O . Note that (| 13[) does not yet imply that the 
every elementary imset is indeed an effective hyperplane defining IC(N), or 
equivalently every elementary imset is an extreme ray of K.*(N). Note that 
an extreme ray of a polyhedral convex cone is a half line and when we say 
"u is an extreme ray of a cone" , we actually mean that {cu \ c > 0} is an 
extreme ray of the cone. At this point we establish the following lemma. 

Lemma 2.1. Each elementary imset is an extreme ray of the imset cone 
IC*(N). 

There are many ways to prove this lemma. We give a somewhat involved 
argument, which will be used frequently in Section |6l 

Proof. Let u = uu^ic) S be an elementary imset. Suppose that u 

is written as a non- negative combination of elements of K,* (N) . By (fl3l) 

«<o,6|Cf>= k v >0, Vv e£(N). (14) 

ves(N) 

We need to show that k v = for all v ^ u^ i, | c) ■ Suppose that for k v > 
for some v — U( a ij,/\c') u (a,b\c)- Among such w's, choose w max = 
u (a'.b'\c) sucn that a'b'C is the largest according to the graded reverse 
lexicographic order of V(N). If a'b'C > abC, then on the left-hand side of 
(fT4|) U( a ,b| c)(a'b'C) = 0, but on the right-hand side of ([T3]l J2ve£(N) k> ' 
v(a'b'C') > 0. This is a contradiction. Hence if k v > 0, then a'b'C < abC. 
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Similarly, considering the graded reverse lexicographic order of C, we have 
C > C, if k v > for v — U( a , b , | c >y In particular for v — U( a / b , | with 
fc„>0we have |C| < \C'\, \a'b'C'\ < \abC\. But this implies '\C\ = \C'\. 
Then for any S such that |5| = |C| we have J2veE(N) ^» ' — 0- Since 
u (o,6|c)('S0 = for S 7^ C, we have fc^ = for i> = u/o'^'is), S ^ C, 
|5| = |C|. Similarly jfc„ = for « = u (a ,, 6 , | C), aW + abC, \C'\ = |C|. 
Therefore > only if C = C and a'b'C — abC. However this implies 

v = u (a,b\C)- □ 

Let C*(N) C R 2 ' 1 denote the linear subspace spanned by the semi- 
elementary imsets. C*(N) is the orthogonal complement of C(N) and its 
dimension is 

dim£*(iV) = 2^1 - \N\ - 1. 

IC*(N) C C*(N) and the dimension of the relative interior of }C*(N) is the 
same as dim£*(iV). By (gj, 

ue£*(iV) ^ u(J5) = 0, ^ u{E) = 0, Ve e JV. 

EeV(N) eeEeV(N) 

Now let 

£J = Z{« <AiB | C> | (AB|C) eT(/V)} 

= { X! fc (A,B|C} • U {A.B\C) I fc (A,B|C) G Z } C Z 2 ' 

(A,B|C)eT(iV) 

denote the submodule of Z 2 ' 1 generated by the semi-elementary imsets. 
Again by (|lip. £J is generated by the elementary imsets: 

C* z = ZS(N) = { K-v\k v eZ}. (15) 

veS(N) 

£J coincides with the set of integer points in C*(N): 

£* z = £*(N)nZ 2lNl . (16) 

This can be seen as follows. For each C C iV, |C| < \N\ — 2, choose an 
elementary imset u/ 0l b|c>- Consider a submatrix of the configuration con- 
sisting these elementary imsets, which is of size 2^1 X (2^1 — \N\ — 1). If C's 
are ordered by the graded reverse lexicographic order as in Table [TJ then 
the lower (2^1 - \N\ - 1) x (2^1 - |JV| - 1) part of the matrix is an upper 
triangular matrix with 1 on the diagonal, as seen in Table [T] Therefore it 
is unimodular (= has integral inversion). Then the projection of C% onto 
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these 2l Jv l-|JV|-l coordinates coincides with Z 2 I N I " 1 . Then (16]) holds, 
because other elements are determined by orthogonality to C(N). 

Let Q + £ (N) denote the set of non-negative rational combinations of 
elementary imsets. An integer point in Q + £ (N) is called a structural imset. 
Let S(N) denote the set of structural imsets. Since K.(N) and fC(N)* are 
rational polyhedral cones, we have 

S(N) = /c*(iV)nz 2lJV| . 

By COI, S(N) can also be written as S(N) = K.*{N) n £|. 

A non-negative integer combination of elementary imsets is called a 
combinatorial imset and 

C(N) = 1+£{N) = { kv -v | k v G Z+} 

t)G£(JV) 

denotes the set of combinatorial imsets. C(N) is the semigroup generated by 
£{N). Clearly C(JV) C 5(iV) and it is known that for \N\ > 5 the inclusion 
is strict (Hemmecke et alP), i.e. for \N\ > 5 the semigroup C(N) is not 
normal. 

3. Multiinformation and derivation of conditional 
independence using imsets 

Let P, /j, be two probability measures on a sample space X ', such that P is 
absolutely continuous with respect to /i. The relative entropy H(P\fi) of P 
with respect to /i is defined by 

H(P\fi) = [ ln^O) dP(x). 

fl"(P|/i) > and H(P\fi) = if and only if P = fj,. 

Consider a joint probability distribution P of the variables in N. As 
usual, A _LL £? | C [P] denotes the conditional independence statement of 
variables in A and in B given the variables in C under P. The set of 
conditional independences under P is denoted by 

M P = {(A,B\C) G T(iV) | AALB\C[P}}. 

We call p the conditional independence model of P. 

For A C N, P A denotes the marginal distribution of the variables in 
A. Let IlieA denote the product of one-dimensional marginal distri- 
butions in A. Then the multiinformation function mp : V(N) —> [0, oo] is 
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defined by 

m P (S) = H(P s \1[PM), 5^0, (17) 
ies 

and mp (0) = 0. Throughout this paper we only consider P such that mp(S) 
is finite for every S C N. The basic fact on the multiinformation is that 
nip is supermodular. Furthermore the following equivalence between the 
conditional independence A _LL B \ C [P] and the local modularity of mp at 
(A, B | C) holds: 

K«(A,B|C))>0, V(A,B\C) eT(N) 

= & A1LB\C[P]. (18) 

This equivalence is the basis for manipulating conditional independence 
statements in terms of imsets. Note that mp is standardized as a super- 
modular function since mp(S) = for < 1. 

Traditionally, the implications among conditional independence state- 
ments under a probability measure P have been studied in terms of the 
following semi-graphoid axioms. In the axioms A, B,C, D C N are disjoint. 

1. triviality A1L$\ C [P] 

2. symmetry A1LB\C [P] => B 1L A \ C [P] 

3. decomposition A _LL BD | C [P] A _LL B \ C [P] 

4. weak union A1LBD\C[P] A _LL L> | BC [P] 

5. contraction A 1L D | BC [F] and A _U_ B \ C [P] A _1L | C [P] 

Note that decomposition, weak union and contraction can be combined 
into the following single equivalence: 

AALBD\C[P] & A1LD\BC[P] and AALB\C[P}. 

This equivalence can be proved very easily by imsets. Take the inner product 
of (|10j) with the multiinformation mp. Then 

(m P , U( A>BD \ C )) = (mp,u {A ^ B \c)}+ (mp, U( A<D \bc))- 
Since every term is non-negative, we have 

A _LL BD \C[P}^ (m P , u {AiBD , c> ) = 

& (mp,U( AtB \c)) = (mp,U( A , D \ BC )) = 
AALD | BC [P] and A JL B \ C [P\. 

In this way, manipulation of "rules" such as the semi-graphoid axioms is 
translated to linear algebraic operations in terms of imsets. This is a very 
important advantage of the method of imsets and the multiinformation. 
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As another example, consider the following identity. 

U(a,b | c) + U(a,c | d) + U<a,d | 6) = i^abc + S c - S ac - Sbc) + (Sacd + Sd - S a d — S c d) 

+ {Sabd + 5b - S ab ~ S bd ) 
= (Sabc + S a cd + Sabd) + ($b + $c + Sd) 

- (Sab + S ac + S ad ) - (S bc + Sbd + S c d) 
= u (a.c | 6) + u (a,d | c) + u (a,b \ d) ■ (19) 

By the method of imsets we then have 

aALb\c,aALc\d,aALd\b <^ a ALc\b,a ALd\c,a ALb\ d, (20) 

where for simplicity we wrote a _LL b \ c instead of a _1L b \ c [P] . It is evident 
that none of decomposition, weak union and contraction can be applied to 
the left-hand side nor to the right-hand side of this equivalence. Therefore 
this equivalence can not be derived from semi-graphoid axioms. 

Here it is interesting to note that (flU)) can be derived by linear algebraic 
operations from (1101) , which itself corresponds to the semi-graphoid axiom. 
By applying (|10p twice we can write 

3 • U( a: bcd | 0) = (U(a,b \ 0) + U(a,c \ b) + U(a y d \ be) ) 

+ ( u {a,c\<$) + u (a,d\c) + U( ay b\cd)) 

+ i u (aA\$) + u {a,b\d) +U(a,c\bd)) (21) 

as well as 

3 • U( a: bcd | 0) = (U(a,b \ 0) + | b) + U( a ,c \ bd) ) 

+ ( u {a,c\ 0) + u {a,b\ c) + U( ay d\bc)) 

+ i u {a,d\$) + u {a,c\d) + u {a,b\cd))- (22) 

Equality of right-hand sides of (|2~T|) and (|2"2"|) immediately gives (|T5| . 

Let u s S(N) be a structural imset. A conditional independence with 
respect to u is defined as follows. We say that A and B are conditionally 
independent given C with respect to u, and denote it by A _LL B | C [u], if 

3k e N, k ■ u - u {A:B | C ) G S(N). (23) 

The set of conditional independences induced by u is denoted by 

M u = {(A,B\C) e T(N) \ AALB\C[u]}. (24) 

We call M. u the conditional independence model of u. The importance of 
this definition lies in the following completeness theorem for imsets. 
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Theorem 3.1 (Theorem 5.2 of Studeny 1 ). Let P be a probability mea- 
sure over N with finite multiinformation. Then there exists u G S(N) such 
that Mp = M u . 

Proof. Define u = Y^,/a B\c)eM P u {a.b\c) - Then clearly Mp C M u . Now 
for any (A, B\C) G M u choose k such that ([23]) holds. Then 

= (m P ,k • u) = (m P ,k ■ u - U( A ,b\c)) + (mp,u^ B \c))- 
Since two terms on the right-hand side are non-negative, we have 
{nip, u^a.b | C)) — an( i hence M u Q Mp. □ 

We have reproduced this proof from Section 5 of StudenjO} because this 
argument is instructive. In this proof, to show (mp, u^ A p \ c)) = 0, we have 
added an extra term k ■ u — uia,b | c) to U/ A ,B \ C)- I n this sense, this proof 
is similar to the argument concerning ([2T|) and ([22]) above for showing (|20|) . 

Next we show that the conditional independence model M u depends 
only on the face of IC*(N) which contains u as a relatively interior point. 
Although this fact is discussed in a series of papers by Studeny^l and Chap- 
ter 5 of Studeny^, we make this point clear in the following lemma. 

Lemma 3.1. Let a\,...,a n be nonzero vectors in M. m . Let C = 
M + {ai, . . . , a n } be the convex cone generated by these vectors and suppose 
that C is pointed and a±, . . . ,a n are extreme rays of C . For a face F of C , 
let F° denote its relative interior. Put Ip = {i \ a, G F}. Then, b G C 
belongs to F° if and only if the following two conditions hold: 

(1) For every i £ If, there exists ki G N such that ki ■ b — a; G C , 

(2) For every i ^ If and every k G N, k ■ b — a,; C. 

Proof, b G C can be written as b = Y^i=i c * a »' c» > 0, « = 1, . . . , n. Suppose 
that b satisfies the above two conditions. Then for i G If, by condition 1, 
for some fcj G N, hb — 4- Xw=i c) a h — 0> or 

1 - 
b = —a l + yj(5j/fef)Oj. 

ki 

3 = 1 

Taking the average of the right-hand side over i G If, we see that we can 
take Ci > for i G If in b — y\ —1 ciai. On the other hand, for i Ip, 
condition 2 implies k ■ Ci — 1 < for any k G N, hence ci — 0. Therefore if 
b satisfies two conditions then b is written as 

b = 2J c % a i, Ci > 0,i e If- 
iei F 
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Since the coefficients are positive, we have b G F°. 

Converseiy let t G F°. For i G Ip and sufficiently large k G N, let a = 
(l/k)cii. Then we have b — a G F, which is equivalent to k ■ b — aj G F C C. 
Let u G K m be a normal vector of a supporting hyperplane of F such that 
(u, x) > 0, Vx G C, and for 6 G F and any a, g" F the inner product with 
respect to v satisfy 

(v,b)=0, (v,ai)>0. 
Then for any k G N, we have 

(w, k ■ b — a,i) < 0, 

which means that k ■ b — a, S" C '. □ 

Under the assumptions of Lemma I3.lt let b G C and let F = F^ denote 
the face of C such that b G F°. Let 

£ b = {di | fe • b — a,i G C, 3k G N}. 

Then £{, is the set of extreme rays of F^ by Lemma 13.11 

From this fact we can prove that conditional independence structures 
induced by imsets depend only on faces of the imset cone K,*{N) and not 
on each imset. 

Proposition 3.1. Let u,u' G S(N). Then M. u = A4 U > if and only ifu,u' 
belong to the relative interior of the same face of KL* (N) . 

Proof. For v G S(N), let F v denote the face of IC*(N) such that v belongs 
to the relative interior F° of F v and let £ v (JV) denote the set of elementary 
imset U( Q .b | c) such, that 3k G N, k ■ v — U( a ^ \ c) G S(N), that is, 

£ V (N) ={« (Bl 6|o I (a,b\C)eM u }. 

As note above, £ V (N) is the set of extreme rays of F v . Hence 

M u = M u , & £ U (N)=£ U ,(N) 

F U = F U ,. □ 

By Theorem 13.11 and Proposition 13. 1\ each conditional independence 
model M.p corresponds to a face of JC*(N). Note that, by mutual orthog- 
onality, the face poset of IC(N) and that of JC*(N) are isomorphic (when 
inclusion in IC(N) is reversed). Hence each conditional independence model 
Mp also corresponds to a face of K,(N). 
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Also note that we can regard the set of conditional independence mod- 
els Aip as a poset with respect to the inclusion relation M.p C A4pr. 
By Proposition 13.11 this poset is a sub-poset of the face poset of K.*(N). 
Unfortunately it is known that for |iV| > 4, there exist faces of K*{N) 
which do not correspond to any Aip, i.e. for \N\ > 4, the set of conditional 
independence models is a proper sub-poset of the face poset of JC*(N). 

Example 3.1 (Example 4.1 of StudenjJS). For N = abed consider 

u = u {c,d\ab) +U(a,b\<$) +U(a,b\c) + U{a,b | d) 

— S N + 2 • 5 ab - 8 ac - 8 a d - he - Sbd - 5 a - S b + S c + 5 d + 6$. 
Then (c, d \ ab), (a, b \ 0), (a, b \ c), (a, b \ d) S Ai u . Furthermore let 

m = A-5 N + 2- 6 s+ Ss - ( 25 ) 

SCAf : |S|=3 SCN:\S\=2,S^ab 

(see Fig. \^) Then it can be checked that {m,v) > for every v S £(N). 
Hence m is supermodular. Furthermore (m,u) = 0, (m, Uuf, i c( j\) = 1. 
Consequently for any k £ N, (m, k ■ u — u^ a _ \ = —1 and hence k ■ u — 
u (a.b\cd) ^ S{N) and (a,b\cd) ^ M u - On the other hand it can be shown 
that for any probability measure P with finite multiinformation 

{c,d\ab),{a,b\9),(a,b\c),{a,b\d) e Mp => (a,b\cd)eMp 

(cf. Corollary 2.1 of Studeny 1 ). Therefore for the above u, there exists no 
P such that Ai u — Mp. 

It can be shown that m in (|25l) is an extreme ray of the supermodular cone 
IC(N) for | TV | = 4 in the sense of the next section. 

4. Results on extreme rays of the supermodular cone 

In this section we prove that certain supermodular functions are "extreme 
rays" of the supermodular cone JC(N). Since extreme rays of K.(N) are 
hard to describe, even partial results on the extreme rays are of interest. 
As we discussed in Section [U since IC(N) contains the linear space jC(N) 
of modular functions, in order to consider "extreme rays" of fC(N) we need 
to identify two supermodular functions / and g if they differ by a modular 
function. For this purpose it is simplest to consider the cone lCi(N) of 
standardized supermodular functions in (jB]): 

K t {N) = {/ G K{N) | f(S) = for \S\ < 1}, 

which is a pointed polyhedral cone in R 2 ' ' ' ^1 1 . 
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Following the terminology of StudenjC} we call an extreme ray of /Q(iV) 
a skeletal supermodular function. Furthermore by abusing the terminology, 
we call / 6 IC(N) a skeletal supermodular function (or an extreme ray of 
/C(iV)) if its standardization / in ([6]) is a skeletal supermodular function. 
Below we often omit "supermodular" and simply say that / is a skeletal 
function, or / is skeletal. Skeletal functions for \N\ = 4 are depicted in 
the appendix of Studeny et alP^ Various conditions on skeletal functions 
of this section are meant to be useful for understanding these functions. 
Except for the extreme ray m in (|25|) . other extreme rays for \N\ = 4 are 
covered by the results of this section. However they are far from enough for 
understanding most of skeletal functions for \N\ =5 studied in Studeny et 
alP 

The results of this section concern extending a supermodular function 
on V(N) to a larger base set N D N and resemble results in Chapter 2 of 
Topkis 11 , although Topkis^ does not consider skeletal functions. Results 
on a skeletal function as a maximum of a collection of modular functions 
are given in Rosenmiiller and WeidnerP^ 

For proving that a supermodular function / is skeletal we proceed as in 
Lemma 12.11 Express / as a sum of two supermodular functions 

/ = h + /a- (26) 

Here we can assume that /, /i,/2 are standardized, namely, = f{S) = 
fi(S) = f2(S) for \S\ < 1. This follows from the fact that the standardiza- 
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tion ([6]) can be performed to /, f\ , f 2 separately. We need to show that f\ 
and fi are proportional to /, i.e., 

/i = cif, fi = c 2 f, ci,c 2 > 0, ci+c 2 = 1. (27) 

We first state two simple facts on skeletal functions. The first one con- 
cerns complementation of subsets of N. In Section 9.1.2 of Studeny^ this is 
referred to as a reflection. Let A c denote the complement of A C N. Let / 
be a supermodular function. Define / by 

f{S) = f(S c ), SCN. 

Taking the complement corresponds to looking at the configuration Un in 
Tabic Q] upside down. The following lemma is trivial. 

Lemma 4.1. f c is skeletal if and only if f is skeletal. 

The second one concerns a supermodular function f(S) depending only on 
the size \S\ of S. 

Lemma 4.2. Let \N\ > 2 and 1 < k < \N\. f(S) = max(|5| - fc, 0) is 
skeletal. 

Proof. It is easily verified that / is supermodular. Let f — fi + f 2 , where 
/11/2 are standardized and supermodular. Since f\,f 2 are non- negative, 
h = h = holds for \S\ < k. This means ([27) holds for |5| < fc. 
Next we consider the case for \S\ > fc. If fc = \N\ — 1, then (|27[) holds. 
Hence assume fc < \N\ — 2. Let Sk be any subset of N such that \Sk\ = fc- 
Let § D S k . Then f{S) = f{T U S k ) where T = S \ S k . For S k C TV, 
define / Sfc (T) = /(T U S k ) = \T\ , T C AT \ S/.. In the same way, define 
f~4 k (T) = fi(T U Sfe), T C N\ Sk, i = 1,2. Then / 5fc and /f fc are modular 
function on N \ S k . By f 1 k and / 2 fe can be written as 

/f * (t) = Af* + £ Af * , ; 2 Sfe (t) = + E » Sk ■ 

From / Sfc (T) = |T| and / = /1+/2, we have Ag* = /i^ = and Af fc = 
1 for a G T. If Af fc = A^ fc = A Sfc and = ^ fc = /A hold for all 
a,/3 G N \ S k where \ s ",fj, Sk G R+, then /f fc (T) = A 5fc |T| and / 2 Sfc (T) = 
/Lt s ' t |T| hold. Furthermore we obtain the representations fi(T) = X\T\ 

and fi k {T) = (jt\T\ where X,fi G R+, because /f fe ("') = if* (a) holds 
if a'S k = aS' k where S' k is a subset of N such that |5fc| = fc and a G 
Sfc, a' G S^. This means that (|2T[) holds for |S| > fc. It is suffices to show 
that Af fc = A?* holds for all a, /3 G iV \ Sfc such that a/ft because we 
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can prove the equation /if fc = (i? in the same manner. Let a, (3 £ N \ Sk 
such that a ^ /3. For 7 £ S*, let 5£ = {S \ 7} U a and 5[f = {S \ 7} U /3. 
First we consider /f fc . From 75^. = aSfc, we have /j k ( r y) = /f fc (a) = Af fc . 
Furthermore, from /3jS' k = a(3S k , we have ff k (£7) = /f* (a/3) = Af* + A^ fc . 
Since / x fc is modular on AT \ S k , we obtain 

/? 08) = fx'" 087) + /? (0) - /f 1 (7) = Af fe . (28) 

Next we consider ff*. From <yS' k ' = (3S k , we have ff k ' (7) = /f fc (/3) = A^ fc . 
Furthermore, from a^S^f = a(3Sk, we have /f fc (07) = /f fc (a/3) = Af fc +A^ fc . 
Since f 1 k is modular on N\S k , we obtain 

jf' (a) = /f * (07) + /f * (0) - /f * (7) = Af * . (29) 

From (gSD, (2g) and /3S£ = aS£ , we have Af fc = /f*(/3) = /f* (a) = Af fc . 

□ 

For stating other results on skeletal functions, we let N = {1, . . . , n}, 
n = \N\, and consider / as a function of xi,...,x„, where Xi £ {0,1}, 
i = 1, . . . , n. S C N can be identified with (xi, . . . , x n ) £ {0, 1}", where 
Xi = 1 if and only if i £ 5. Let Aj/ denote the function depending on Xj, 
j ^ i, defined as 

\f = f(xi,- ■ ■ ,Xi_i, l,X i+ i, . . . ,X„) - /(Xi, . . . ,Xj_l,0,X i+ l, . . . ,x„). 

Then it is easily shown that / is supermodular if and only if 

A,A,./ 0. Vi^i. 

Note that AjAj-/ is a function of Xfe, fc ^ Furthermore / does not 
depend on Xj if and only if A^/ = 0. Also since a standardized supermodular 
function / is non-decreasing, Aj/ > 0. 

From now on we interchangeably write f(S) or /(xi, . . . , x„), i.e., the 
argument of / may be a set or a 0-1 vector, depending on the context. We 
now prove the following lemma. 

Lemma 4.3. Let A C N and let g : V{A) -> R. Define f : V(N) R by 
f(S) = g(SC\A). Then f is supermodular if and only if g is supermodular. 
Furthermore f is skeletal if and only if g is skeletal. 

Proof. We omit the proof of the first statement, since it is trivial. Also it 
is easily shown that if g is not skeletal, then / is not skeletal. It remains to 
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show that if g is skeletal, then / is skeletal. Write / = /1 + /2, where /, fa, fa 
are standardized. Then = Ajf for j £ A implies = Ajfa = Ajfa, j A. 
Hence /1 , fa do not depend on Xj ,j $ A and /1 , fa can be considered as 
supermodular functions on V(A). Then since g is skeletal, fi,fa have to be 
proportional to g = f. □ 

For AC N let 

/a(5) = Uc(S) - 1 { ac5} = r lf f" 5 (30) 

I otherwise. 

Using above lemmas we prove the following proposition. 
Proposition 4.1. Ja = 1ac« *s skeletal. 

Proof. On £>(A), = g A (S) = max(|S| - |A| + 1,0) is skeletal by 

Lemma 321 Then its extension to V(N) is skeletal by Lemma T4. 31 □ 

We now consider more complicated extensions of supermodular func- 
tions. For / : {0, 1}" K define 

/o(a;i,...,x n _i) = f(xi, . . . ,x„_i, 0), 
fx(xx,. . . ,x„_i) = f(xi, . . . ,x„_i, 1). 

The following lemma is easy to prove. 

Lemma 4.4. Let f be a standardized supermodular function, such that 
fa = 0. Then f is skeletal if and only if fa is skeletal. 

The following result is somewhat more difficult to prove. In the lemma 
and its proof we denote N' = { 1 , . . . , n — 1}. 

Lemma 4.5. Consider f such that fa is a skeletal function satisfying 
A Jo{N' \ i) = 1 for all i E N' and fa(S) = \S\, S C N'. This f is 
skeletal. 



Proof. We first show that / is supermodular. It suffices to show that 
AiAjf > for all i,j € N,i ^ j. For i,j ^ n and S C N' \ ij, 
AiAjf(S) = AiAjfa(S) > because of the supermodularity of fa. Also, 
for i,j ^ n and n € S C N \ ij, A i A i /(S' \ n) = A t Ajfa(S) > 0. Thus 
we need to show the supermodularity when j = n and S C N'. From the 
supermodularity of fa, for i £ N' \ S we have 

fa(Si) - fa(S) < fa(N') - fa(N' \ i) = 1 = fa(Si) - fa(S), 
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which means that 

= A t (f 1 (S) - f Q (S)) = Aifx(S) - AJ (S) > 0. 

Therefore / is supermodular. 

Now write / = g + h where g, h are standardized. Define go,gi,h±,ho 
as above. Considering the — 0, since /o is skeletal, there exist 

c g , Ch > 0, c g + Ch = 1, such that 

.9o = Cgfo, h = c h fo- 

Also from the assumption and supermodularity of g\ and h\ , for any i e N' 
we have 

gi (N') - 91 (N' \ i) > g (N') - g (N' \ i) = c s A 2 / (iV' \ i) = c g , (31) 
fti(iV') - hi(N' \ i) > h (N') - h Q (N' \ i) = c h A t f (N' \ i) = c h . (32) 

If at least one of the inequalities is strict, we have 

1 = (n - 1) - (n - 2) = h(N') - h{N' \ i) 
= (gi(N') - 9l (N' \ *)) + (hi(N') - h!(N' \ *)) 
> c g Aif (N' \ i) + c h Aif (N' \ i) = 1, 



which is a contradiction. Therefore ([31]) and ([32]) are equalities. 

Now consider the case that x„ = 1. fi(S) = \S\ is modular. Hence both 
gi and h\ are modular as well. From ([31]) and the modularity of g\ , for any 
i G TV' we have 

c 9 = gi(N') - . 9l (TV' \ i) - <&(») - 5l (0) = 

Furthermore from /i(0) = and the modularity of <7i, we have 51 (0) = 0. 
Hence by (|3J), we obtain 

91 (S) = J>(<) - c 9 |S| - cJ^S), 5 C TV'. 
Therefore g = c 9 /. Similarly h = ctf- This proves that / is skeletal. □ 

For our next result, we write £c„_i = {x%, . . . , a;„_i). Let /' : {0, 1}" — > 
M be supermodular. Define / : {0, l} n+1 -> M by 



/ («^n— 15 f) 1? 

/' (a;„_ i,0) otherwise. 



Then the following result holds. 

Proposition 4.2. J/ /' is skeletal, then f is skeletal. 
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Proof. As in the last proposition we write N' = {1, . . . ,n — 1}. We also 
write N" = {1, ...,n + 1}. We first show that / is supermodular, i.e. 
AiAjf > for all i, j E N", i ^ j. By assumption, A t Ajf > for i,j E N'. 
Next consider the case i 6 N' and j = n + 1. Then by definition 

A n+ i/(x n _i, 1) = A n f'(x n -i), (33) 
An+i/fx,,-!, 0) = /'(asn-i, 0) - /'(aj n -i, 0) - 0. (34) 

Hence 

A i A n+1 /( i r n _ 1)5 , 1) - AiA n /'(x n _ l | ) > 0, AiA„ + i/(a3 n _ li j, 0) = 0, 

where x n l j = aj^yw. The case of i G A^',j = n is similar. For i — n, j = 
n + 1, by taking the difference of (|34|) and (p3|) we have 

A n A„ + i/(x n _0 = A n /'(x n _i) > 0. 

Therefore / is supermodular. 

Now let / = g + h. Repeating x n define 

f {x \ , . . . , X n ) — f \X\ , . . . , x n , x n ) . 

Then / = /'. By assumption / is skeletal. Hence defining g, h similarly, we 
have 

g = c g f, h = c h f, c g , c h > 0, Cg + c h = 1. 

This implies that g(x) = c g f(x) and h(x) — Chf(x) if x n — x n+ \. Therefore 
we only need to consider the cases (x n ,x n+ i) = (1, 0) or (x n , x n +i) = (0, 1). 
By symmetry we consider only the former case. Then 

c g /(cc n _i,0) =g(x n -i,0) = g(cc„_i,0, 0) < g{x n -x, 1, 0). 

The last inequality holds by monotonicity of g. Similar relation holds for 
h. If at least one of the inequalities for g and for h is strict, then 

/(x„_i, 0) = c g f(x n -i, 0) + c h f(x n -i, 0) 

< g(a;„_i, 1,0) + /i(se„_i, 1,0) = /(aJn-1, 1,0) 
= / / (x„_ 1 ,0)=/(x„_i,0) J 

which is a contradiction. Hence 

c g f(x n -x, 1,0) = c ff /'(x n _i,0) = c s /(a; n _i,0) = g>(x n _i,l,0) 

and g = c g /. Similarly = Cft/. Hence / is skeletal. □ 
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Our final result of this section concerns a product of skeletal functions 
for disjoint sets. These functions played an important role in Kashimura et 
aP 

Proposition 4.3. Let AB = N, A n B = and consider f of the form 

f(S)=g(AnS)h(BnS) (35) 

where g : V{A) — > K and h : V(B) M. Suppose that f,g, h are standard- 
ized supermodular functions. Then f is skeletal if and only if g and h are 
skeletal. 

Note that if g and h are standardized supermodular functions on V(A) 
and V(B), respectively, then f(S) in (|35p is a standardized supermodular 
function on V(N). 

Proof. It is easily shown that if g or h are not skeletal, then / is not 
skeletal. Hence if / is skeletal, then both g and h are skeletal. 

We now show the converse. Suppose that g and h are skeletal. Write 
/ = f\ + f 2 , where /i, f 2 are standardized. Consider any A' C A such that 
g(A') > 0. For S C B, consider 

f(A'S) = MA'S) + f 2 (A'S) = g(A')h(S), SCB. 

Since fi(A'S), f2(A'S) as functions of S are supermodular and since h is 
skeletal, there exist ci(A'), C2{A') > 0, such that 

h(A'S) = ci(A')h(S), h(A'S) = c 2 (A')h(S). 

Hence 

f(A'S) = a{A')h{S) + c 2 (A')h(S) = g(A')h(S), S c B. 

Note that this holds also for A' such that g(A') = by defining = 
ci(A') = C2{A'). Now fixing any S C B such that h(S) > and considering 
f(A'S)/h(S) = g(A') = Ci(A') + c 2 (A') as a function of A', we see that ci 
and C2 have to be proportional to g because g is skeletal. Then f\ and f 2 
are proportional to / = g/i. □ 

In Appendix A, we present a generalization of Proposition 4.3 for arbi- 
trary cones. 
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5. Structure of faces of semi-elementary imsets 

In this section we study a face of the imset cone corresponding to a semi- 
elementary imset uia,b | c) ■ By the face corresponding to Uia,b | c) > we mean 
the unique face of the imset cone, such that Uia,b \ c) is m the relative 
interior of the face. We denote this face by F( A B | ^ . In Kashimura et alP 
we have shown some remarkable facts on F^ A B | c y Here we establish more 
basic facts on F^ AB | c ^ . 

If we repeatedly apply the decomposition in (1101) , any semi-elementary 
imset M(a,b|c> can De written as a sum of elementary imsets. Those ele- 
mentary imsets are from the following set of elementary imsets: 

£(a,b\c) ={u{a,b\v) \ ae A,beB,C crc ABC}. 

In this section we establish the following basic facts on F( AB \ c) '■ i) £(a,b \ c) 
is the set of extreme rays of F( A ,b | C); u ) dim(F( AB \c)) = (2' A ' — l)(2l B l — 
1). For proving these facts we give linearly independent set of supermodular 
functions orthogonal to Fi A ,b \ c) ■ 

Let D — N \ {ABC), i.e. N = ABCD. For E C N define l E c. by ([30]) . 
which is skeletal by Proposition 14. II Similarly define 



^ EDS 
otherwise, 



which is skeletal by Lemma 14.11 since 1ed» = ljc c ,' Now consider the 
following supermodular functions. 

Uxcc. {A l C A), (36) 

l Sl cc. (0 ^ C B), (37) 

1bc i2 . (£? C AB, Ci c CO, (38) 

1bd x c. C ABC, ± D 1 C £>). (39) 

Let Arabic) denote the set of these supermodular functions. The cardi- 
nality of M(a,b | c) is given by 

\M (AtB \c)\ = 2 |A| + (2 |s| - 1) + 2l AB l(2l c l - 1) + 2 l ABC l(2l D l - 1) 

= 2 |JV| _( 2 pA| _ 1 ^ 2 \°\ -1). (40) 

First we check that the linear independence of the above supermodular 
functions. 



Lemma 5.1. The elements of M{a,b\c) are linearly independent. 
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Proof. We need to show that the coefficients ,^b ± , ^ec\ , £e.Di £ R in 

^ ^A 1 'i-A 1 CC»(S) + ^2 MBi lfii CC.(S') + ^ ^£Cil£CiD.(S') 

AiCA 0^BiCB BCAB.CiCC 

+ X! £WbdiC.(S) = (VS) 

are zeros. Let S — A 2 C (A 2 C A). Then Xm iC ^ 2 = ^ anc ^ by induction 
we have = 0. Similarly by letting S = B 2 C (0 ^ B 2 C _B) we have 
S0#b 1 cb 2 MBi = and by induction = 0. 
Next let 5 = £ 2 C 2 (E 2 C AB, C 2 cC). Then 

XI = 

E 2 C£CAB,C 2 CCiCC 

and the double induction on E and Ci yields KeCi = 0- Finally letting 
S = E 2 D 2 (E 2 C ABC,0 D 2 C D) we have Ebcb^b^ 6^ = 
and the double induction on E and Z?i yields C-E_Di =0. □ 

Now we prove the following theorem. 
Theorem 5.1. An elementary imset u belongs to £(a.b\c) if an ^ only if 

(«,/) = 0, VfeM (AiB \o. (41) 

Proof. We first show the orthogonality (u, f) — for all u G £(a,b \ c) an d 
for all / G M {a ,b\c)- Consider f(S) = l Al cc. (Ai C A). Let 

u (aM \C>) = <5a6C" + fc' - - foe £ £<A,B | C) • 

Then 

<«<a,6|C')./> = fW) + f(C) - f(aC') - f(bC') 

= ^-{A 1 CCabC} + IjAiCCC"} ~ ^-{A 1 CCaC} ~ ^-{A^'CbC} 
= 1 {A 1 CCaC'} + ^-{A 1 CCC} ~ ^-{A 1 CCaC} ~ ^{A 1 COC} 
= 0. 

Similarly f(S) = 1b 1 cc, is orthogonal to u< ,6|c") € £<a,b|c>- 
Next consider f(S) = l E dD. (E Q AB, C\cC). Since C % C x 

(u {aM c>),f) = f(abC) + /(C) - /(aC) - /(&C) 

= ^{ECiDabC"} + 1{£CiDC'} _ ^{ECi^aC"} ~ ^-{ECiDbC} 

= + 0-0-0 = 0. 
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Finally consider f(S) = 1 E d 1 c. (E C ABC, ^ D x C D). Then 

<«<»,& 1 CO , /) = f(abC) + /(C) - /(aC) ~ f(bC) 

= l{BDiCa6C"} + 1{£DiCC'} — !{F£>iCaC"} — i {ED 1 CbC'} 

= + 0-0-0 = 0. 

We have shown that (u, f) = 0, Vu G £<a,b i c>> V/ G .M/y^s | c) ■ 

We now show that for any elementary imset u = u/ a / ^ | c") ^ £{a,b | c) 
there exists / G ■M(a,b\c) sucn that (it, /) = 1. For / = Ifc. we here 
confirm which elementary imset U( a ,^|r) satisfies (u( Ql( g | r) , /} = 1- Write 

ri = rn-F, r 2 = rnf c . 

Then we can ignore T2 and replace T by Ti in 

lFc Q/ 3r + Ifcr — lFc»r — 1-Fc^r- (42) 

Hence we only consider Ti. It is easy to see that (|42p is equal to 1 if and 
only if the first term of (|42|) is 1 and other three terms are zeros. In this 
case 

h8ri|>|F|, |r x |<|F|-2 

implies F — a/3Ti. Putting T2 back we have the following equivalence 

(«(a,|S|r>,lFc.} = l a,P€F, (af3T) n F = F 

& a, 13 G F, af3T D F. (43) 

Similarly for 1fd» we can show 

(u(a,p I r> , Ifd.) = 1 ^ a.jE'F, TEF. (44) 

Now arbitrarily fix a, f3, T. Consider the case a(3T % ABC , i.e. (a/3T) n 
D ^ 0. Let 

D ± = (aPT) n D, E = (a(3T) n ABC. 

Then ED X = a/3T and a,/3 E F, a/3T D F, where F = ED X . Hence 
(|43|) implies (u/ aj p \ r), IfDiC.) = 1- It remains to consider the cases where 
a/3T C ABC. Hence assume C % T C ABC. Let d = CnT. Then Ci C C. 
Let B = m (AS). Then T = Ed. For F = T, we have a, /3 £ F, T = F. 
Hence (j44|) implies (it( QiJ a|r)i 1fCiD») = 1- It remains to consider the case 
r D C, a/3r C ABC. Note that we want to eliminate the case that both 
a and /3 belong to A or belong to B. By symmetry we can only consider 
the former case. Hence assume a, f3 E A, T D C. In this case we can set 
A x =aj3(Anr). Then ((43]) implies (u^ir), UiCC.) = 1. □ 
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Remark 5.1. In the above proof we have shown that 

(w,1ac.),(«,Ud.) G {0,1} 
for every elementary imset u and every A C N . 

Remark 5.2. Let Li^A^^s | c)) and Lin(f^ B | C .)) denote linear sub- 
spaces spanned by M.ia,b \ c) an d £{a.b \ C)i respectively Linearly extend- 
ing (|41l) we see that Lm(A4 i A ,b \ c>) is the orthogonal complement of 
Lin(f( j4 B | C )). 

Two facts stated at the beginning of this section are simple consequences 
of Theorem 15.11 

Corollary 5.1. £(a.b\c) * s ^ e se ^ of extreme rays of F(a,b\c) an d 
dim(^ (A , B | C) ) = (2l A l-l)(2l B l-l). 

The first statement is obvious and the second statement follows from 
Remark 15.21 and (|40|). The fact that &voh(F(a,b\c)) does not depend on 
C can also be seen from the one-to-one linear correspondence between 
Lm(£(A,B | c)) and lim(£( A>B \ >) given by 

£{A,B\C) 3 u {a,b\T) <r+ U( a< b\r\C) G £(A,B\$)- 

Note that Theorem 15.11 shows that elements of M.i A ,b\C) are extreme 
rays of the face Fh b\c) °^ tne supermodular cone K(N) dual to F^ AB | c) ■ 
However it is hard to believe that F^ b\c) ^ s siniplicial and hence F^ A g | c ^ 
should have other extreme rays than elements of M.ia,b \ c)- Characteriza- 
tion of extreme rays of Fh^ b \q^ is an interesting question left to our future 
research. 

Remark 5.3. In our previous manuscript Kashimura et alP we studied 
how the semi-elementary imset U(a,b \ c) is expressed as a non-negative 
integer combination of elements of £(a,b\c)- By Theorem 15.11 if Uia,b\C) 
is expressed as a non-negative integer combination of all elementary imsets, 
then the coefficients of u ^ £(a,b \ c) have to be zero. Therefore Theorem 
I5.1l iustifies our restriction to imsets from £(a,b | c) f°r expressing uia,b \ c)- 

6. Small relations among imsets 

Consider the configuration Un as in Table [1] The kernel of Un is denoted 
by kerWjv. Furthermore let 

ker z Z^v = (ker U N ) n l) £ ^ 
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denote the integer kernel oHAn- Let z = {z u ) u e8(N) € be an element 

of the kernel of Un- Considering / in ([5]) as a row vector we have 

= fU N z = (l,l,...,l)z= 2 z*. (45) 

Hence any non-zero z G keri/jv has both positive and negative elements. 
We denote the positive elements of z as at, . . . , and negative elements 
of z as — j3i, . . . , —/?m- Then (|45|) can be written as 

aiui H h cik u k — Pi v i H h^mUm, (46) 

where ui, . . . , u/s, vi, . . . , v m are distinct elementary imsets. We call this 
kind of equality between two non-negative combinations of elementary im- 
sets a relation among elementary imsets. More precisely we call (j4"6)) a k 
by m relation. Since — z also belongs to kerWjv, we can take k < m. Also 
we can assume that u\, . . . , Uk and v\, . . . , v m are ordered according to the 
graded reverse lexicographic order of elementary imsets. We say that the 
relation (|4*S| contains another relation J2i'=i a 'i' u 'i' = S?=i P'j' v 'j' ^ ev ~ 
ery u' it appears on the left-hand side of (|46[) or every v'j, appears on the 
right-hand side of P5|) . 

The most basic relation is the following 2 by 2 relation, which comes 
from the semi-graphoid axiom. 

U(a.M | C) + u (a.b 2 biC) = u (a,b 2 \ C) + u (aM \ b 2 C) ■ (47) 

The corresponding element of ker Un is written as 

8{a,bx | C) + &(a,b 2 | bxC) ~ S(a,b 2 | C) ~ S(aM \ b 2 C) i (48) 

where the unit vector <5( a ,fc|c) S l) £lyN ^ is defined as <5( a ,6 1 c) ( w ) = 1 if 
v = u (a,b\c) an d otherwise. In connection to Markov bases discussed in 
Section [7] we call an element of ker Wat a move. Let B denote the collection 
of moves in (|48p over all a, 61,62, C. We first show that B generates the 
integer kernel of Un, i.e. B contains a lattice basis for Un- 

Theorem 6.1. ker^l^N as a submodule ofZ> £ ( N '> is generated by B. 

The following proof also shows that B spans keri^v as a linear subspace 
ofRl^WI. 

Proof. Let AT = {1, ... , n}. Let IB be the submodule of Z' £ ( w ) generated 
by B. Since the inclusion ker^Wjv 2 TLB holds, we prove kev-^UN Q ZS. 
Consider the graded reverse lexicographic order of £ (N) in ([9]). Recall that 



27 



we always assume a < b if we write uu t b \ c) an d the maximum element of 
£(N) with respect to < is u/ 1)2 1 ^ , where R = {3, . . . , n}. 

For any non-zero z € Zl £ WI W e define the degree deg(z) by the mini- 
mum elementary imset ^ | q\ such that 2(tf/ ^ c) ) 7^ 0- Now fix any non- 
zero z e keiz^N- We will find some v G ZZ5 such that deg(z — v) > deg(z) 
or z — v. Then the proposition follows by induction. We first show that 
deg(z) is not the maximum u« 2 im. Indeed, if deg(z) = uh2\r\, then z 
has to be z(u(i^\r))&{i,2\r) an d therefore 

= U N ■ z = z(u<i, 2 |fl))u( li2 |fl), 

which contradicts to z{uh2\r)) ^ 0. Hence deg(z) < uh^IR)- Denote 
deg(z) — M( a .b|c) an d let £c — {uu'^no) I a ' \b' S N\C}. Denote the 
maximum of £c by U( a ,p \ c)- We show that U( a j, | c) < u {a,/3 \ c)- Indeed, if 
u {a,b | c) — u (a.p | c) i then z can be written as 

z = z{U(a,p\C))${a,P\C) + X! X! z ( u (a',b' \ C'))^(a',b' \ C) ■ 

OC a',b'£N\C 

Then 

= U N ■ Z = z(U( ai p\ C ))U( aj p\C) + ^2 X! Z ( U (a'M'\C')) u (a'M'\C) 

C>C a',b'€N\C 

= z (u {a ^\ c) )8 c + (terms of {S C '}c>c), 

which contradicts to z(v,( a p | c)) 0- Hence U/a^ic) < u/ a ^ic)- 

Now we give a vector v e Z£> such that deg(z — v) > deg(z) or z = v. If 
b < (3, then consider a vector 

V = ^(W( ,6|C)) {8{a,b\C) + $(a,/3 | bC) ~ &{a,f3 \ C) ~ ^{a,b\/3C)) 

in ZB. Note that deg(w) = W( a ,b|c) = deg(z) and the leading terms of v 
and z are the same. Therefore we have deg(z — v) > deg(z) or z — v = 0. If 
b = (3, then a < a. Consider a vector 

V = z{u( a ^ b \ C )) (5{ a ,b\C) + 5{a,b\aC) ~ 3{a,b\C) ~ &(a,b\aC)) 

in "LB. Note that deg(u) = W( a ,6|c) = deg(z) and the leading terms of v 
and z are the same. Therefore we have deg(z — v) > deg(z) or z — v = 0. 
This completes the proof. □ 

For the rest of this section we consider k — 2, 3 in (J46J) and show that 
only certain types of relations exist. By (|4"5j) we have 



ax H h ak — Pi H h /3 m 



(49) 
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Our first result on small relations is the following. 

Theorem 6.2. Consider a relation in (|46|) with k — 2: 

Ot\U\ + a 2 U 2 — fav\ H 1- PmVm, (50) 

where oti, a 2 , fa, ■ ■ ■ , /3 m are integers. Then m — 2 and oti = o>2 = fa = fa 
and the relation is a positive multiple of (147[) . 

Proof. Let u = a x ui + a 2 u 2 . From (j49|) . we have 

a x + a 2 = fa + ■ ■ ■ + fan, ax > 0, a 2 > 0, fa > 0, . . . , fa n > 0. 

Because u±,u 2 ,vi, . . . ,v m are extreme rays of JC*(N), there is no relation 
like axUi + a 2 u 2 — favi- Therefore, we have m > 2. Furthermore, be- 
cause ui, U2, vi, . . . , v m are distinct elementary imset, the following proper- 
ties hold: 

ui^u 2l Ui,u 2 y^Vi,...,v m , Vi^vj (51) 

for i ^ j. We write the elementary imsets in (|50[) as Ui = U( ai bl Cl ),it 2 = 
u (a 2 M I c 2 ) , Vi = "(ai.bi |c;).---,fm = u (a' m .b' m \C' m )- Remember that the 
terms are ordered such that u± < U2 and v± < V2 < ■ ■ ■ < v m according 
to the graded reverse lexicographic order of elementary imsets. Then there 
exists li G {1, . . . , m} such that C\ = C[ = ■ ■ ■ = C[ and Ci +1 , . . . , C' m ^ 

d. Note that if C\ ^ C 2 , then a x = fa H 1-/3^ and di&iCi ^ a-fr-C- for 

i G {1, . . . , h}. In the same way, there exists I2 S {1, . . . , m} such that 

fl2 6 2 C 2 = a' h b' l2 C"i 2 (52) 

andC 2 ^C/ 2 . 

Now we define the weight w for a relation as follows: 

w = 2_j u(S) = 51 (aiUi(S) + a2U 2 (S)). 

{SCN I |S| = |Ci| + l} {SCN I |S| = |Ci| + l} 

Then the possible values of w are — 2ai, — 2«i — 2a 2 and — 2«i + a 2 . 

Firstly, we consider the case of w = —2ai. Then we have |C 2 | > \Ci\ + 2. 
In this case, ui — Vj holds for some j € {1, . . . , m} and this contradicts (|5Tj) . 
Therefore there is no relation with w — —2ai. 

Next, we consider the case of w = —2cei — 2a 2 . We prove Ci ^ C 2 . If 
Ci = C 2 , then by (l52|) there exists some ^ S {1, . . . , m} such that vi = u 2 . 
This contradicts (f5Tj) . Hence we obtain Ci 7^ C 2 . In the same way, we can 
prove aibiCi ^ a 2 6 2 C 2 . Now, since \C X \ = |C 2 | in this case, we have C 2 = 
C[ , j = • • • = C T ' n . However, there exists Z 2 e {1, . . . , m} such that ([52]) 
holds. If I2 € {?i + 1, . . . , m}, then this contradicts ([51]) . If I2 £ {1, • • • , h}, 



29 



then C\ C a 2 b 2 C 2 and C2 C a\b\C\. This means that C\,C 2 C aibiCi PI 
a 2 b 2 C 2 . Then we have Ci U C2 C a^Ci n a 2 b 2 C 2 . From Ci 7^ C2 and 
|Ci = C2I, we have \C\ UC2I > \C\\ + 2. Furthermore, from aibiCi ^ 
a 2 b 2 C 2 , we have \a\b\C\ f~l 0262^^2 1 < \aib\Ci\ — 2 = |Ci|. However, this 
means that |Ci| +2 < |Ci|. This is a contradiction. 

Finally, we consider the case of w = —2u\ + a 2 . Let ft — \C±\. Then we 
have |C 2 | =h + l, |oi6iCi| = ft + 2, |a 2 6 2 C 2 | = ft + 3 and |Cf| = ft + 1 
for i e {lx + 1, . . . , to}. This means that a 2 b 2 C 2 = a ; 1 +i^; 1 +iC;' 1 +i = ■ • • = 
a' m b' m C' m . Now, assume u(axbiCx) > 1. Since Vj(5) < for 5 = ft + 2 and 
ie{!i + l,..., to}, we have ai&iCi = a'^bjCj for some j G {1, . . . , 1%}. This 
contradicts (|5T|) . Therefore u{a\b\Ci) < and we have 

> u(a 1 biCi) = axU\{a\bxCx) + a 2 u 2 (aibiCi) = oc\ + a 2 u 2 (aibiCi). 

From the above equation, we have u 2 (a\biCi) < — ^ < 0. This means 
u 2 (aibiC\) = — 1. Therefore, we have a 2 > u\. In the similar way, we 
obtain u(C 2 ) < and a.\ > a 2 . Hence we have ct\ = a 2 . From this, we have 
u(C 2 ) = u{aib\Ci) — 0. Then there exist d±, d 2 S {ai, bi} and c?3 € {02, 62} 
such that 

C 2 = Ci U {di}, O161C1 = C 2 U {d 2 } = Ci U {di,d 2 }, 
a 2 6 2 C* 2 = aihd U {d 3 } = C* 2 U {d 2 ,d 3 } = C x U {di.da.da}. 

Then, from u(CiU{d 2 }) = w(a 2 6 2 C 2 \{e? 2 }) = —1, there exist i' 6 {1, . . . , l{\ 
and j' S {1% + 1, . . . , m} such that 

v v {Ci U {d 2 }) = v f (a 2 b 2 C 2 \ {d 2 }) = -1. (53) 

Let i>i' (Ci U CZ4) = —1 for some g?4 6 N\{CiUd 2 }. Then there exist j' such 
that Vji{C\ U <&i) = 1. Since Vj'(a 2 b 2 C 2 ) = 1 for this j', we have C\ U c/4 C 
a 2 b 2 C 2 = C% U {di, g?2, e^}. This means 1^4 G {^1,^3}. Furthermore, since 
7^ Wj' for i' 7^ j' in (|5"Tj) . tty is limited to one of Vi> = M<d 2j d 3 1 d) or 
v i' = u (d 3 ,d 2 1 Ci) j according to d 2 < c/3 or ^3 < <i 2 . In the same way, Vf 
is limited to one of U( dl ,d 2 1 d 3 d) or M(d 2 , c j 1 1 d 3 Ci), according to dt < d 2 or 
d 2 < d%. Therefore we obtain m = 2. Furthermore from (f53|) . the relation 
is a positive multiple of (j47|) . □ 

One consequence of this theorem is the following corollary on two- 
dimensional and three-dimensional faces of K* (N) . 

Corollary 6.1. Two-dimensional faces of K.*{N) are simplicial cone gen- 
erated by two arbitrary elementary imsets which do not appear on the same 
side of the relation (|47j) . Three-dimensional faces are either simplicial cones 
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having three elementary imsets as extreme rays or a cone with four extreme 
rays appearing in the relation (1471) . 

Proof. Two-dimensional cones are always simplicial. If two elementary im- 
sets do not span a two-dimensional face, then the two-dimensional cone 
spanned by these two imsets cuts a relative interior of a face of a higher 
dimension. Then it intersects a cone generated by other extreme rays of 
this face. Hence these two imsets appear on one side of a 2 by m relation. 
Then by Theorem 16.21 to — 2 and this face has to be of the form (|47|) . If a 
three-dimensional cone is not simplicial, then any set of four extreme rays 
are linearly dependent. Hence there is a relation among them, which has to 
be a two by two relation. By the same argument it is easily seen that there 
are no more than four extreme rays of a three-dimensional face. Then by 
Theorem 16.21 this face has to be of the form (|47|) . □ 

Another important example of relation is (|19l) with k = to = 3. 

Theorem 6.3. Consider a relation in (|46|) with k — 3: 

aiui + a 2 u 2 + a 3 u 3 = PiVi -\ h/3 m w m , (54) 

where oi\, a 2 , a 3 , j3\, . . . , (3 m are integers. Then one of the following two 
properties holds. 

(1) to = 3, ot\ = a 2 = «3 = Pi — fi 2 — @3 and the relation (fl>4")) is a 
positive multiple of 

u (a,b 1 | b 2 C) + u (a,b 2 | b 3 C) + u (a,b 3 | biC) 

= u (a,b 2 I fciC} + u (a,b 3 I b 2 C) + u (aM | b 3 C) (55) 

(2) The relation (|54[) contains a relation of Theorem \6.lX 

The statement in (1) can not be strengthened because there exist 3 by 
777 (to > 3) relations of type (2) such as 

2U(a,c | 0) + U(a,b | c> + U(b.c \ a) = ^(a.b | 0} + u (b,c \ 0) + ^ u (a.c \ b) 

(= U(ab,c\0) + u (q,6c|0)), 

and 

u (a,b 3 | fciba) + 2M( a ,fc 2 | + u (aM \ 0) 

= U(a,b 2 \bib 3 ) + u (a,b 3 \bi) + u (aM \ b 2 ) + U{a.b 2 | 0)- 
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Proof. Let u = oi\U\ + a 2 u 2 + ct 3 u 3 . From (|49j). we have 

(Xi + ot2 + «3 = /3i H 1- fei ai > 0, a 2 > 0, o: 3 > 0, /?i > 0, . . . , /3 m > 0. 

Because Hi, 1/2,^3, i>i, ... , v m are extreme rays of }C*{N), we have to > 3. 
Furthermore, because Ui, 112, v\, . . . , u m are distinct elementary imsets, the 
following properties hold: 

ui^u 2 , ui,u 2 ^v 1 ,...,v m , Viy^Vj. (56) 

We write the elementary imsets in ([M]) as iti = u^ ai bl | , u 2 = 

U {a 2 M I C 2 >) u 3 = "(a 3 ,b 3 I C 3 )! U 1 = U {a' 1 ,b' 1 \ C[) , ■ ■ ■ : v rn = u (a' m M' m C' m ) ■ As in 

the proof of Theorem 16.21 we define the weight w for a relation as follows: 
w = 2J (aiui(5) + a 2 u 2 {S) + a 3 u 3 (S)). 

{SCN I |S| = |Ci| + l} 

Remember that the terms are ordered such that u\ < u 2 < u 3 and 
Vi < v 2 < ■ ■ ■ < v m according to the graded reverse lexicographic order 
of elementary imsets. Then the possible patterns are classified as follows: 



(i) 


w = 


-2(ai +a 2 + a 3 ) <*\Ci\ = 


= IC2I - 


\Ca\, 


(ii) 


w — 


-2(ai +a 2 ) + a 3 \Ci\ = 


= IC2I - 


|C 3 [ — 1, 


(iii) 


w — 


— 2ai + a 2 + a 3 \C\ \ — 


\C2\-l- 


= \c 3 \ - 


(iv) 


w = 


-2«i + a 2 \Ci \ = \C 2 \ - 


- 1 < \c 3 


-2. 


(v) 


w = 


-2ai & \Ci\ < \C 2 \ - 2 < 


\C3\-2. 




(vi) 


w = 


-2(a 1 + a 2 )^|Ci| = |C 2 ] 


< \c 3 \ - 


- 2 



First, we consider the case (jvj) of w — —2u\. Because u\ ^ Vi from (52), 
there exists i' G {1, . . . to} such that = (ai, d! \ Ci), d! ^ bi. To cancel out 
the value iv (d'Ci) = — < at d'Ci, there exists i" e {1, . . . to} such that 
Vi» = (a" ,b" I d'Ci), d' / 61. Furthermore there exists d" e {a", b"} such 
that Ui(d'd"Ci) = for i = 1,2,3. To cancel out the value Vi„(d'd"Ci) = 
—oti" < at d'Ci, there exists iv» = (a"',b"' \ d'd"C\). In the same way, 
there exists d!" S {a'", b'"} such that Ui{d' d" d'" d) = for i = 1,2,3. 
Since it is impossible to cancel out the value Vi>»(d l d"d l "Ci) = —ai>>> < 
at d'd"d"'Ci, this case can be ignored. 

The case (jvi) is almost the same as in the case of (jvj). 

We consider the case (jrvj) . If | C2 1 — 1 < | C 3 1 — 2 , then it contradicts (j56|- 
Therefore, we consider the case of \C 2 \ — 1 = \C 3 \ — 2. In this case, there 
exist mi and m 2 such that 

l<TOi<m 2 <TO, ai=Pi-\ +/3 mi) Ci = C' 1 = --- = C' mi , 

"3 = Ati 2 +i H 1" Pm, a 3 b 3 C 3 = a' m2+1 b' m2+1 C' m2+1 = ■■■ = a' m b' m C' m , 
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and we have \C^\ = \C m + i\ = • • • = \C m2 1. Now let a\b\Ci = {di,da}UCi. 
From ([56]) and the same argument of the case (jvj), we can write C*2 = d\C\ 
without loss of generality. Furthermore from (|56[) . we have a\ < ai. If 
u(a\b\C\) < 0, then aiiti + aiu-i has a relation of Theorem l6.2l Therefore, 
we have u{a\b\C\) > 0. In the same way, we have u{C 3 ) > 0. Hence there 
exist i' G {m2 + 1, . . . , m} and j' G {1, . . . , mi} such that C-/ = ai6iC*i and 
a'jb'jC'j G (7 3 . This means aihCi ^ C 3 . Since C\ C did = C 2 C ai&iCi = 

C a-,6-,(7-, = 03636*3, we can write 03636*3 — d 2 d 3 d^C2 = d^d^d-sd^Cx. 
From ([56]) . we have diC\ C C3. This means that C3 is equal to did 3 C\ or 
did±C\. Without loss of generality, we assume that C3 = d 2 d 3 Ci. In the 
same way, we can show that 02626*2 is equal to 03636*3 \ d\ or 03636*3 \ 04. 
Since 026262 D C2 = d\C\, we have 026262 — 0363(73 \ c/4 = did2d 3 C\ 
and ^(diG^Ci) = — 1. However, from aifeiCi = di^Ci and ai < 02, 
zi(ai6iCi) = ai — a2 < 0. This is a contradiction. 

Next, we consider the case dm}. If u(C2), u{C 3 ) > 0, then a contradiction 
follows from the proof of Theorem 16.21 Therefore, either u(6*2) < or 
u{C 3 ) < holds. Without loss of generality, we assume that ^(62) < 0. 
From ([56]) . we have u(ai6i6i) < 0. Therefore, we obtain 



> u(aib 1 C 1 ) = ai + a2U2(ai6iCi) + a 3 M3(ai6i6i) 

OL2 Ot 3 

^=^> — u 2 (aib 1 Ci) H u 3 (ai6iCi) < — 1. 



From this we have 1*2(016161) = —1 or ^3(016161) = — 1. However, if we 
assume ii2(oi6i6i) = — 1, then a\U\ + a 2 U2 has a relation of Theorem 16.21 
Therefore, we have 112(016161) — and 113(016161) = — 1. From this, we 
have C2 C 0161C1 C 0363(73. In this case, if 1/(63) > 0, then a contradiction 
follows from the proof of Theorem l6.2l Therefore 1/(63) < holds. Further- 
more, we have ui(Cs) = — 1. This means that otxiii + 0:31/3 has a relation 
of Theorem O 

The case ([n]) is almost the same as in the case of (fm|) . 

Finally, we consider the case Since (7- G {C\, C2, C3}, a^6-(7 4 ' G 
{016161,026262,036363} and (f56|) . we have m = 3 and (7- = (7j. In this 
case, a'^o' x C' x = 02626*2 or Oi6'i6( = 03636*3 holds. Without loss of gen- 
erality, we assume a' 1 b / 1 C[ — 0262(72. Then, we have a^b^C^ — 0363(73, 
03636*3 = ai6iC*i and u\ = ct2 = a 3 — j3\ = 02 = fa. Therefore the 
relation is a positive multiple of the three by three relation in (1551) . □ 
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7. Some computational results on Markov basis for imsets 

In this section we give some computational results on Markov basis for the 
configuration Un- Before giving them we discuss motivation for studying 
Markov basis from the viewpoint of inferences on conditional independence 
statements. In Theorem l6.1l we showed that the set of moves in (|48|) contains 
a lattice basis for ker^ Un- This means that if we allow negative coefficients, 
then all relations with integer coefficients can be derived from the two by 
two relations of Theorem 16.21 On the other hand, recall that the three by 
three relation in Theorem 16.31 can not be derived from two by two relations 
in the sense of our discussion after (|20|). This is because in "applying a 
rule" , we do not allow subtracting a non-existing term from either side of 
a relation. If a Markov basis is available, then we can derive any relation 
without ever subtracting a non-existing term starting from a given initial 
relation. For definitions and details of Markov basis refer to Chapter 5 of 
Sturmfels^or Chapter 1 of Drton, Sturmfels and SullivantP^ 

Let Q(Un) be a (minimal) Markov basis for Un and \Q(Un)\ be the 
number of elements in Q{Un)- If \N\ < 2, then there is no relation in Un 
and hence kerWjv = 0. If |iV| = 3, then there is only the relation of (|47|) with 
C = 0. We now give some computational results on Markov basis for N = 
4, 5. In the following results, we identify two relations if one is obtained by 
permuting characters in another. For example, we identify the following two 
relations: U( a<b \ d ) + W( 0)C | bd) = U( a>c \ d ) + U( a ,6 1 cd) and U( a<b \ c ) + u< M | ac) = 
u {b,d\c) + u (a,b\cd)- Let G r (Un) be the set of representative elements in 
Q{Un) under this equivalence relation. Let |^y fi (Wjv)| be the number of 
elements in Q r {Un)- 

Let M(a,b | c) be a sub-configuration of Un, N = ABC, which consists of 
all of the elementary imsets u^jj \ p) such that a £ A, 6g_B, CCTC ABC, 
a < b and rna = rnfe = 0. Markov basis for Uia,b \ c) represents the re- 
lations which appear in the face F^ B | c ^ of IC*(N). Markov bases for the 
subconfiguration of elementary imsets belonging to F^_ B | c ^ is of some in- 
dependent interest because they generate a combinatorial pure subring in 
the sense of Ohsugi, Herzog and Hibf^H^. Let Q R ((A, B | C)) be the repre- 
sentatives of Markov basis for U^.b \ c) an d |^? fl ((^4, B\C))\ be the number 
of elements in Q R ((A,B\C)). Let Q R (T(N)) = J2abc=n S r ((A, B \ C)) 
and |£? fl (T(iV))| be the number of elements in G R (T(N)). 

For a relation of (|46|) . we call a\ + • • • + ajc the degree. Note that 
Ei=i a i = E"LiPrn- We give g R {U N ) and G R {T(N)) at each degree for 
| N | = 4, 5 in Table [2] and El The Markov basis of Un were computed with 
4ti223. Since the computation has not completed when we used the graded 
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reverse lexicographic order in ([9]), we give the result for the following order: 

h c >2 c \ 

u (a,b | c) < "(a',fe' | c) if < or C = C and b > b' , 

I or C = C and b — b' and a > a', 

where 2 C = £ j£C 2\ 

Table 2. Number of representatives of elements of Markov basis for |JV| = 4 

degree g R (U N ) Q R (T(N)) 

2 2 2 

3 1 1 

4 4 4 

Table 3. Number of representatives of elements of Markov basis for \N\ = 5 
degree G R (U N ) Q R (T(N)) 



2 


3 


3 


3 


2 


2 


1 


11 


11 


5 


18 


16 


6 


210 


162 


7 


384 


36 


8 


364 


38 


9 


90 





10 


220 





11 


16 





12 


63 






In Table [2] for \N\ — 4, all of the representatives in G r (Un) appear in 
G R (T{N)). However, in Table[3]for \N\ — 5, there are many representatives 
which appear in G r (Un) but do not appear in G R (T(N)). Especially there 
is no relation with degree more than 8 in G R (T(N)). This means that the 
facial structure of IC*(N) is quite complicated. Furthermore, through this 
computational study, though there are a lot of coefficients more than 1 in 
the Markov basis for G R (T(N)) for \N\ = 5, we confirmed that the Markov 
basis for G R ((A, B \C)),\A\ + \B\ + \C\ < 5 is square free i.e., the coefficients 
in the relations are all 1. It is future work to confirm whether this is true 
or not for |iV| > 6. 
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Appendix A. 

Here we present a generalization of Proposition 14.31 to convex cones on a 
direct product of two finite sets. 



Proposition 7.1. Let X,Y be finite sets and let G C {g : X — > 

H C {h : Y — > R + } be pointed convex cones in M x and R y respectively. 
Define 

M(G, H):={f:XxY^R+\ Vx, f(x, •) e H, Vy, f(; y) e G}. 

Then M(G,H) is a convex cone. If g and h are extreme rays of G and H, 
then g{x)h{y) is an extreme ray of M{G,H). 

Proof. Let g and h be extreme rays of G and H and write g(x)h(y) = 
fi(x,y) + f2(x,y) (/i,/2 £ M(G, H)). Note that there exist x' , y' such that 
g(x') > and h(y') > 0. For x' with g(x') > we have 

,/ x fi{x',y) }2{x',y) 

Kv) = — rrr + — rrr> Vy - 

g(x') g{x') 

Since h is an extreme ray of H and the two terms on the right-hand side 
belong to H as functions of y by the definition of M(G,H), fi(x',y) = 
ci(x')h(y), f2(x',y) — C2(x')h(y) for some c\(x'), C2(x'). For x with g(x) = 
0, we have fi(x,y) — fi{x,y) — (Vy) and we can put c\(x) = C2(x) = 0. 
Then for all x,y we can write fi(x,y) = c\(x)h(y), f2(x, y) — C2(x)h(y). 
In particular for y' with h(y') > 0, c\(x) = fi(x,y')/h(y'), C2(x) = 
fi{x, y')/h(y') belong to G as functions of x and 

fl (,) = ^4g^ = c l(a; ) + C2 (,). 

Hence c\ 1 C2 are proportional to g. □ 



We now show that Proposition l4.3l is a special case of Proposition[7jTJ We 
note that skeletal functions are also extreme rays of K,+ (N) := JC(N) D {/ 
/ > 0}. In fact, if / e )C e (N) and f = fi + f 2 (fufi 6 K+(N)), then it is 
easily seen that /i,/2 £ Ki{N). 

Let X = 2 A ,Y = 2 s , G = fC + (A), H = K.+ (B), and let g and h be 
skeletal on A and B. Then g and /i are extreme rays of G and H . Hence by 
Proposition 17.11 f is an extreme ray of M(G,H). Now from the definition 
of supcrmodularity K t {AB) C fC+(AB) C M{G,H). Therefore / is an 
extreme ray of ICe(AB). 
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Remark 7.1. We can not put G = lCe(A), H = Ki{B) in the above proof, 
because Kt{AB) % M(G,H). For example, let f(S) = l{ a fccs} (a E A,b e 
B), then / G K t {AB) but f(AS) = l{bes} (VS C B) docs not belong to 
and/£M(/Q(A),/Q(-B)). 

Acknowledgment. We are grateful to Milan Studeny for very detailed 
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